Dynamical generation of gauge groups 
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It has been known that the dynamics of k coincident D-branes in string theory is described effectively 
by U(fc) Yang-Mills theory at low energy. While these configurations appear as classical solutions in 
matrix models, it was not clear whether it is possible to realize the k 1 case as the true vacuum. 
The massive Yang-Mills-Chern-Simons matrix model has classical solutions corresponding to all the 
representations of the SU(2) algebra, and provides an opportunity to address the above issue on a 
firm ground. We investigate the phase structure of the model, and find in particular that there exists 
a parameter region where O(N) copies of the spin-f/2 representation appear as the true vacuum, 
thus realizing a nontrivial gauge group dynamically. Such configurations are analogous to the ones 
that are interpreted in the BMN matrix model as coinciding transverse 5-branes in M-theory. 
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Introduction. — The discovery of D-branes as classi- 
cal solutions in string theory suggested an interesting sce- 
nario that we are actually living on a kind of D-brane em- 
bedded in a higher-dimensional space-time ("the brane- 
world scenario" ) . The low energy effective theory of a 
D-brane is given by a U(l) gauge theory, and in the case 
of k coincident D-branes, the gauge group enhances to 
U(fc). Along this line, there are a lot of activities in 
the search of a perturbatively stable brane configuration 
which realizes the Standard Model at low energy. 

On the other hand, it is also possible that our world 
is realized as the true nonperturbative vacuum of super- 
string/M theories. Such an issue may be addressed by 
studying the matrix models [1,2], which are proposed as 
nonperturbative formulations of superstring/M theories. 
In particular the IIB matrix model [2] can be obtained 
by dimensionally reducing the 10-dimensional M — 1 su- 
per Yang-Mills theory to a point. The space-time should 
come out dynamically as the four extended directions in 
the eigenvalue distribution of the 10 bosonic matrices [3]. 
The mechanism that realizes the brane world in this way 
has been studied in Refs. [4,5], and by now there are cer- 
tain evidences that indeed /cwr-dimensional space-time is 
generated dynamically [6-9]. 

The possibility of obtaining the U(fc) gauge group dy- 
namically has been discussed in Ref. [10]. There it was 
claimed that if the eigenvalues of the ten bosonic matri- 
ces form clusters of size k, the low energy effective theory 
should have the U(fc) gauge symmetry. However, whether 
such clustering really occurs as a dynamical property of 
the IIB matrix model remains unclear. 

In the series of papers [11-13], we discussed the same 
issue in matrix models with a cubic term. In these mod- 
els the fuzzy spheres [14] appear as classical solutions 



[15] and play the role of D-branes [16]. One of the ad- 
vantages of studying fuzzy spheres in the present con- 
text is that they are solutions even for finite matrices, 
and therefore one can study their dynamical properties 
by well-defined perturbative calculations. (See Refs. [17] 
for other motivations for fuzzy spheres.) The expansion 
around the k coincident fuzzy spheres yields a U(fc) gauge 
theory on a noncommutative geometry [15]. Therefore, 
by comparing the free energy for various fuzzy sphere 
configurations, we may discuss the dynamical genera- 
tion of gauge groups. However, in the simplest Yang- 
Mills-Chern-Simons matrix model [11,12], as well as in 
its higher-dimensional extensions [13], which accommo- 
date four-dimensional fuzzy manifolds, the gauge group 
generated dynamically turned out to be U(l). 

In this letter we show that in fact one can obtain gauge 
groups of higher rank dynamically by adding a "mass 
term" [18] to the matrix models. A similar model is 
known to appear in M-theory on a plane wave back- 
ground [19]. Although we consider that such a phe- 
nomenon occurs in more general models, here we discuss 
the "massive" Yang-Mills-Chern-Simons matrix model 
for simplicity. 

The model and its classical solutions. — The model 
we consider in this letter is defined by the action 
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S = N ahi i -- [Ap, A v f + - i e^xApAvAx 



^P 2 {A,) 2 \ , (1) 



where A^ (/i 



1,2,3) are N x N Hermitian matrices. 
(Rewriting the action in terms of A^ — a A^, one obtains 
the model of Ref. [11] plus a mass term \ M 2 (Aft) 2 with 
M = a p.) The classical equation of motion 



1 



-[A v , [Ap, A„}} + i e^x [A v , A x ] + p 2 A l _ l = 



(2) 



has a class of perturbatively stable solutions given by 
A^ = = x L n f OT < p < ^75, where (p = 
1,2,3) is an arbitrary A-dimensional representation of 
the SU(2) algebra and x = § (l + V 1 - 2 P 2 ) ■ Usin S 
the SXJ(N) symmetry, the representation matrix can 
be brought into the block-diagonal form with fcj blocks of 

(n) 

the rii-dimcnsional irreducible representation L\\ satis- 
fying n\ < n 2 < ■ ■ ■ < n s and J2t=i n * = N. (The 
trivial solution A^ = may be regarded as a particu- 
lar case in which one chooses N copies of the represen 



should therefore restrict the matrix C, as well as h, to be 
within the tangent space of H. Since A(X^) is given by 
the product of the eigenvalues of the above linear trans- 
formation in the restricted space, we obtain 



a (n,+n J )/2-l 
i,j=l l=\ni~rij\/2 



x 2 i{i + i) 



k i k j (2l+l) 



(6) 



where the symbol Y[ ' implies that I = is excluded. 
Inserting the unity (5) and exploiting the SU(iV) invari- 
ance, we rewrite the partition function as 



tation L 



(i) 



0.) Since the Casimir operator for each Z = vol(H) / dA8{i[X^ A^] - C) A(A tl ) > 



(7) 



irreducible block is given by (l^^J = \ (n 2 — 1) 1„, 
the configuration X^ may be viewed as a collection of fcj 
coincident fuzzy spheres with the radius \ \\J i n i) 2 — 1; 
and the classical action can be evaluated as 

S c i = \Na 4 f( X ) Y^hm^mf-l) , (3) 



where f{\) = \ X 4 - f X 3 + 5 P 2 X 2 - 

Since f(x) is positive (negative) for p > | (p < |), 
the minimum of 5 c i is given by the single fuzzy sphere 
A li =x. for p < I , and by the trivial solution A^ = 
for p > |. These configurations describe the true vacuum 
in the large-a limit where quantum fluctuations are sup- 
pressed. As we decrease a, however, the one-loop effects 
become non-negligible, and it is possible that some other 
solution describes the true vacuum. 

One-loop calculation. — Let us evaluate the parti- 
tion function Z = J dAe~ s around a general solution 
A^ = Xfi at the one-loop level. The measure for the 



path integral is defined by dA = 11)1=1 Yia=i dA®, where 

A^ = J^ali A l ta with ta bcin g thc generators of \J(N) 
normalized by tr (t a t b ) = 5 a b- 

For solutions other than A^ = 0, we need to fix the 
gauge since there are flat directions corresponding to the 
transformation A^ 1— > A 9 ^ = g A^ , where g is an el- 
ement of the coset space H = U(N)/ U.UA V ( k i)]- We 
take the gauge fixing condition 



i[X li ,A li ] = C , 



(4) 



where C is a Hermitian matrix, and consider an identity 

1= / dg5{i[X„A*\-C)b{AJ , (5) 

where the Faddeev-Popov determinant A(A M ) needs to 
be evaluated for A^ = X^ at thc one-loop level. Making 
an infinitesimal transformation g = 1 + i h on A^ = X^ 
in (4), we obtain [A M , [X M , h]] — C. The linear transfor- 
mation defined by the left hand side has zero eigenvalues 
for h belonging to thc Lie algebra of rii=iU(fcj)- We 



where vol(iJ) = vol(U(iV))/ r]- = i vol(U(/c i )) can be ob- 

tained by vol(U(n)) = 7~j^ k.. q i • Since the result does 
not depend on C, we integrate over it within the tan- 
gent space of H with the Gaussian weight Ne~z Na tr c , 

where Af — ( J . This yields an ex- 

tra term S g .f. = — \ N a 4 tr [X^, A^] 2 in the action, which 
lifts thc flat directions as desired. 

Now we are ready to perform the integration over A^. 
Decomposing the variables as A^ = X^ + A^, we ex- 
pand the total action Stot = S+ Sg.f . with respect to the 
fluctuation A^ as 



S to t = Sci + \ Na 4 tr (ApQ^Aj) 



+ ■ ■ 



where the operator <2 M „ is given by 

= \x 2 {£\) 2 + P 2 } -ip 2 e^\ £\ 



(8) 



(9) 



We have introduced an operator C\ which acts onaJVx 
N matrix M as C\ M = [L\ , M] . Thc remaining task is 
to solve the eigenvalue problem Q^A^ = A A^. Since 
has a block-diagonal form, we can decompose the matrix 
A^ into blocks of size ni x nj, and it suffices to solve 
the eigenvalue problem within each block [20]. This can 
be done in a similar way as in the BMN matrix model 
[21]. As a complete basis for each block, we choose the 
eigenstates \l,m) of (£\) 2 and £ 3 with the eigenvalues 
1(1 + 1) and m, where |m| < / and \rii — nj\/2 < I < 
(n, + nj)/2 — 1. Expanding each block of A^ as 



(A^m 



(n l +n J )/2-l l 

E E 

l=\rii— rij |/2 m=—i 



I, m) 



and using the properties 



(£1 ± 1C2) 1 1, m) = b± 1 1, m ± 1) , 

b± = y/l(l + l)-m(m±l) 



(10) 



(11) 
(12) 



we obtain thc eigenvalue equation within each block as 
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{A-^(l-m)}4' m) =P^-4' 
{A - p 2 (l + m)}A { b m) = -p 2 b+A9' m+ V 



7 (l,m— 1) 



where we have introduced A± m ^ = A^' m ^ ± i A^ and 
A = A — x 2 I {I + !)• For Z > 1/2 we have eigenvalues 



(13) 

3 (14) 

(*,m+l) ; (15) 



Ai(0=X 2 i(* + l) , 
\ 2 {i) = x 2 i{i + i)-P 2 i , 



(16) 
(17) 
(18) 



whose degeneracy is 2Z + 1, 2Z — 1, 2Z + 3, respectively. 
(For 1 = 1/2 the 2nd eigenvalue does not appear.) For 
I = we have A = p 2 with 3-fold degeneracy. Thus the 
integration over A^ yields 
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(n+m)/2-l 



IL 



-3(fc,) 2 



ri(^) M 



'■.;=i 



9n m = I]' [A 1 (/) 2i+1 A 2 (/) 2i - 1 A 3 (/) 

i=|n-m|/2 

Bringing all the factors together, we get 
Z = vo\(H)AfA(X f ,)Ze- s " 1 . 
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(19) 
(20) 

(21) 



In fact this result holds also for the case — formally. 

If we take the large-A r limit of the free energy F = 
— log Z with fixed a and p, the one-loop terms give 
maximally O (N 2 log N) contributions. The 0(7V 2 logiV) 
terms coming from vo\(H ) and Af cancel each other in 
general. Apart from the universal term ^N 2 \ogN com- 
ing from the first factor in (19), we have contributions 
from A(A p ) and the third factor in (19), which induce 
extra 0(Af 2 logA^) terms when m becomes of O(N). 

The p < | regime. — Since the single fuzzy sphere has 
a negative classical action of 0(a 4 A^ 4 ) in this regime, no 
other solutions can be the true vacuum unless a becomes 
as small as 0(^=). The free energy for the single fuzzy 

sphere A^—x is obtained as 



N 2 



FS 



N 2 a A f{ X ) + -log^ + 4 log a -J (22) 



at large N, where 5 is an 0(1) constant which depends 
only on p. (The one-loop calculation is reliable for the 
single fuzzy sphere as far as a^/N is large [11,12]. This 
is not the case for all the solutions.) 

On the other hand, at small a we find, after rescal- 
ing A,j, = aA^, that the quartic term in the action (1) 
becomes dominant, and obtain the so-called Yang-Mills 
phase [11], which is described by the vacuum of the pure 
Yang-Mills model [22]. Using the result of the Gaussian 
expansion method [23], we obtain the free energy of this 
phase as 



J^ f ym = - logTV + 3 loga + 7 



(23) 



at large N, where 7 ~ —4.5. 

Comparing (22) and (23), we find that Fym becomes 
smaller than Fps at 



a < a c 



1 /2 1ogA/V /4 

Tf\~]J(x)\ 



(24) 



This is consistent with the Monte Carlo simulations at 
p = [11], where a first order phase transition has been 
observed with the upper and lower critical points oi^ ~ 
0.66 and ~ respectively. Notice the inequality 

«cr < a cr (p = 0) < oi^- 

The p > I regime. — Since f(x) > in this regime, 
in order for some nontrivial solution to have smaller free 
energy than the A^ = solution, it should have a clas- 
sical action of 0(A^ 2 ) or of smaller order. Thus it turns 
out to be sufficient to consider the case where s and all 
the rii are of 0(1), while ki are of O(N). Such a configu- 
ration (except A^ = 0) is analogous to the ones that are 
interpreted [24] in the BMN model as a collection of co- 
incident transverse 5-branes in M-theory. Let us denote 
hi = Ti N, where are real parameters of 0(1). In what 
follows we consider the large-A limit taken in this way. 

Let us first consider the case where a is very large, so 
that we can neglect all the a-independent terms in the 
free energy and obtain the relevant terms 



jL F „ ni {{n t f l} + B^n) 2 , 



(25) 



where A — \a A f(x) and B — 2 log a. We determine the 
r» that minimize the free energy within the constraint 
J2t=i n i r i — 1 using the Lagrangian multiplier A as 
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A rii — A n 



{M 2 -i} 



(26) 



where A should be fixed by the constraint. We may 
assume n p — p, and requiring > 0, we obtain 
s(s 2 - l)(4s 2 - 1) < GOB /A, where s should be taken 
to be the largest possible integer. 

Thus we obtain the A^ = solution in the large-a limit 
as expected, but as we go below the critical a determined 
by ^ = |, O(N) numbers of 2 x 2 blocks start to appear. 
If we decrease a further, larger and larger blocks appear 
with the specific proportion (26). This result is valid 
when p is very close to | so that the coefficient f(x) in 
A is small and the transitions take place at large a. 

At moderate a, the one-loop terms which are indepen- 
dent of a can no longer be neglected. We search for the 
minimum of the free energy numerically restricting our- 
selves to the case n p = p (p = 1, • • • , s) with s = 7. In 
Fig. 1 we plot the result for p = 0.7. (The situation for 
other values of p(> I) is qualitatively the same.) Since 
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we obtain r-? = throughout the whole region of a, the 
above restriction is considered to be harmless. The min- 
imum free energy obtained in this way is plotted in Fig. 
2, where we also plot the result (23) for the Yang-Mills 
phase, which should be valid for small a, and expected 
to be smoothly connected to the = solution at 
a <~ 0(1). In the p > | regime, since the action (1) 
is positive definite, we can easily prove that the free en- 
ergy should be a monotonously increasing (continuous) 
function of a. We therefore consider that higher loop 
corrections increase the free energy for the "5-brane" at 
a ~ 0(1), where the Yang-Mills phase should take over. 

Combining all the results obtained above, we arrive at 
the phase diagram depicted in Fig. 3. 
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FIG. 1. The parameters (2 < i < 6) that minimize the 
free energy are plotted against a for p — 0.7. 
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FIG. 2. The "free energy density" limjv_» 00 (-^F — | log N) 
is plotted against logo for p — 0.7. The solid line represents 
the result obtained by the one-loop calculation, which is valid 
at large a. The dip corresponds to some "5-brane" solution, 
and the straight part corresponds to the T M — solution. 
The dashed line represents the result (23) obtained for the 
Yang-Mills phase, which is valid at small a. 
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FIG. 3. A schematic view of the phase diagram of the mas- 
sive Yang- Mills- C hern- Simons model. The critical point a CT 
between two phases below p — | is given by eq. (24). 
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